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ABSTRACT 

The partition function of the 2d Ising model with random nearest neighbor cou- 
pling is expressed in the dual lattice made of square plaquettes. The dual model is 
solved in the the mean field and in different types of Bethe-Peierls approximations, 
using the replica method. 
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INTRODUCTION 



The application of methods of mean field type to Ising models allows one to 
obtain very accurate approximations of the thermodynamic quantities. However, in 
presence of quenched disorder this approach is of difficult implementation. In this 
paper, we show that rather good results can be obtained after performing a duality 
transformation of the Ising model with random nearest neighbor coupling that assume 
the values J^- = ±1 with equal probability. The model is thus defined on a dual 
lattice where the spin variables are attached to the square plaquettes. The advantage 
is that the quadratic term of the dual Hamiltonian has constant coefiicients instead 
of random ones. It is therefore possible to use the standard methods of the mean 
field to estimate the quenched free energy. Our results can be generalized to higher 
dimensions, although the approximations become more rough, because the number of 
plaquette spins over the number of interaction links increases with the dimensionality 
[1] . In particular we obtain an extremely good estimate of the ground state energy of 
the random Ising model, by applying the Bethe-Peierls approximation where part of 
the short range order is taken into account. 

In section 1, we introduce the dual lattice made of elementary square plaquettes. 
On this lattice the partition function can be expressed as a function of the inverse 
temperature (3 = — | lntanh(l/T) where T = is the temperature of the original 
lattice. 

In section 2, we apply the mean field approximation to the dual model, using the 
replica method. 

In section 3, we introduce the Bethe-Peierls approximation. This allows us to 
obtain a very precise estimate of the ground state energy of the two dimensional Ising 
model with random coupling. 
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In section 4, we show that it is possible to improve the Bethe-Peierls approximation 
by introducing an interaction between different replicas. 

In section 5, the reader will find some remarks and conclusions. 



1. DUALITY TRANSFORMATION 

The partition function of the Ising models on a lattice of N sites with nearest 
neighbor couplings Jij which are independent identically distributed random variables, 
in absence of external magnetic field, is 

ZnW, {Jij}) = n exp(/?Jy o-jCTj) (1) 
{^} (id) 

where the sum runs over the 2^ spin configurations {a}, and the product over the 2 A'" 
nearest neighbor sites One is interested in computing the quenched free energy 

1 



/ = - lim In Zn (2) 

where A indicates the average of an observable A over the distribution of the random 
coupling. The quenched free energy is a self-averaging quantity, i.e. it is obtained in 
the thermodynamic limit for almost all realizations of disorder [2]. 

On the other hand, it is trivial to compute the so-called annealed free energy 

corresponding to the free energy of a system where the random coupling are not 
quenched but can thermalize with a relaxation time comparable to that one of the 
spin variables. In our model, where the coupling are independent dichotomic random 
variables Jy = ±1 with equal probability, one has 

fa = -r' ln(2cosh2/5) (4) 
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However, fa is a very poor approximation of the quenched free energy, and is not able 
to capture the quahtative features of the model. 

In order to estimate (1), it is convenient to use the link variable Xij = ai aj, since 
only terms corresponding to products of the variables Xij on close loops survive after 
summing over the spin configurations: on every close loop of the lattice Yl^ij — 1; 
while Yi ^ij = CTa C6 for a path from the site a to the site b. A moment of reflection 
shows that it is sufficient to fix J| Xij = 1 on the elementary square plaquettes V to 
automatically fix it on all the close loops. The partition function thus becomes 



where the number of plaquettes is 



2 

{Xij} 1=1 



Np = N 



and we have introduced the plaquette variable Xi = Ylvi 

For dichotomic random coupling = ±1 with equal probability, the free energy 
of the model is invariant under the gauge transformation one 
has 

where Jj = Ylvi again a dichotomic random variable (the 'frustration' [3] of the 

plaquette Vi). It is worth remarking that (6) gives the partition function in terms of a 
sum over the 2^^ configurations of the independent random variables Xij = ±1 with 
probability 

= 2^^ 

In this section we shall indicate the average of an observable A over such a normalized 
weight by 



e.g. one has {xij) = tanh/5 and (xi) = tanh^/?. With such a notation, the partition 
function assumes the compact form 



N 

Zn = 2^ cosh2^(/3) ( J](l + xj,) ) (8) 

In order to estimate the average in (8), let us introduce the dual lattice [4] where 
the sites are located at the centers of each square of the original lattice. A dual spin 
variable is attached to each square plaquette and can assume only the values ai — ±1 
with equal probability, so that one has the identity 

1 + xJi = J2 {^iJii^^~"'^'^- (9) 

Since there is a one-to-one correspondence between links on the original and on the 
dual lattice, we can estimate the link-average noting that 

(E^^ ' ) = (E ' ' ) = E (tanh/3)^ (10) 

The last equality in (10) follows from the fact that 

, (i+g,) < Xij >= tanh/9 if ai 

1 a (Ti = (Tj 



Inserting (10) and (9) into (8) one has 



N 



= 2^ cosh^^(/?) Yl E J^'^"'^^^ E (12) 

{a} i'i'j) 

where we have introduced the variable 

P=--lntanh(3 (13) 

which is the inverse temperature of the dual model vanishing as e~^^ when the 
temperature T = — > 0. The quenched free energy (2) thus becomes 

= lnsinh(2/?) - pf0) (14) 
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where / is the free energy of the dual model, defined as 

/(^) = - hm J-lnZj, (15) 

N^oo 

in terms of the partition function 

From (16) the Hamiltonian of the dual model can be defined via the relation 

^iVp = E{~}e"^-^, as 

iy = -5]?,?, - J]ln(J,)lii^ (17) 

^=1 

Let us stress that the quadratic term of the dual Hamiltonian is independent of the 
random coupling and the randomness enters via a random complex magnetic field 
that can assume the two values and i7r/(2/3) with equal probability. In fact, the 
weight exp(— H.) does not define a standard Gibbs probability measure on the dual 
lattice: it defines a signed probability measure, differing from that one of the pure 
Ising model only for the presence of the random sign related to the frustrations of the 
square plaquettes {Ji}. 

2. REPLICA TRICK AND MEAN FIELD APPROXIMATION 

The introduction of the dual model allows one to apply the mean field approximation, 
since one can easily linearize the Hamiltonian (17) by neglecting fluctuations. A 
similar method has been introduced in the framework of field theory in statistical 
systems without disorder such as lattice gauge theories or spin models [5-6]. 

For our purposes, it is convenient to use the replica method in order to get the 
quenched free energy of the dual model as 

/(^) = - hm hm J- In (Z^ (18) 
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Let us thus consider n non-interacting replicas of our disordered system labelled by 
a = 1, • • • , n. Now, the J's are are independent random variables in 2d. Indeed, one 



can easily verify that Yli Ji = Hi Ji because = ±1 with equal probability. It is worth 
noting that this is not true in Sd, where the Ji of the i square plaquette of a cube can 
be obtained as a product of the remaining five J's of the cube, implying nctt6e Jk = ^i 



so that ricwbe <^fc = 1 while ( Jfc )^ = 0. On the contrary, a plaquette frustration Jj 
cannot be expressed as a product of the other ones in 2d. As a consequence, from (16) 
the partition function of n replicas becomes 



N n 



{s} i=l a=l 

where the sum in (19a) runs over the 2^"^ spin configurations {s} of the replicas, and 
we use the compact notation: 

{s}^{5W},.--, 
One can easily perform the disorder average in (19a) and gets 

= J^/E„E,,,^.^"^^r ;Q 1 (l + i-irHa^A (196) 

{s} i=l \ a J 

As the free energy is invariant under the gauge transformation a^"-* — > —al"\ (19b) 
assumes the simpler form 

jz^=j2'^E,-E,,ri-''r' fi 1 (i+nsjA (19c) 

{s} i=l \ a / 

It is worth stressing that the above expression differs from the partition function of 
a collection of n non-interacting Ising systems without disorder only because of the 
factor Ylii^ + Y[a^i°'^)/'^- Such a term introduces an 'effective' interaction between 
replicas: a configuration contributes to the annealed partition function only if 
ria^i'^^ = 1 on each site of the dual lattice (plaquette of the original lattice). 
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Now we can use the mean field approximation to estimate (19), by introducing the 
magnetizations 



rua = hm o-J"^ a=l,---,n (20) 



Indeed, if we neglect the fluctuations, the quadratic term of (19b) can be estimated 
as CTj^"^ CT^"^ = so that (16) becomes 

{s} i=l \ a / 

The mean field solution can be found by the introduction of n auxiliary fieldS 
$1, • • Using the saddle point method, one has, in the limit A?" — > oo. 



r exp (at 2 /? (2 mA - ^l)) (22) 

J — oo 



As a consequence, (Z^)^ is given by the maximum over ■■■,$„} of 

{s} i=l V a / 



{s} i 



n _ n ~ ^ 

;r(«)<T, -TT —Crv^ ^flrh ^(«) 



||e4/3ar*«+|| 



a 



(23) 



Now, we can explicitly carry out the sum over the 2^"^ spin configurations in (23) and 
obtain 

N ^ / n 

e L^a = \ « 

$1, 



_ N ^ / n n \ 

(Z;^= max e-'^'^SLi^^H ||2cosh(4^$«) + J]2sinh(4^$«) (24) 



In 2d, it is commonly believed that there is no glass transition and no replica symmetry 
breaking, so that we expect that the maximum of (24) is realized at the same $a = $* 
for all the replicas. As a consequence, using the replica trick (18), the quenched free 
energy in the mean field approximation reads 

f0) = max ^ ln(2 sinh(8 /3$) ) - 2/3 (24) 



where the maximum of (24) is reahzed by the value solution of the self-consistency 
equation 



The graphical solution of this implicit equation is showed in fig 1. One sees that 
should be always larger than unity and at /3 — > cxo (infinite temperature T — (3~^ limit) 

= 1. It can appear rather odd that in the dual model the magnetization > 1. 
This stems from the fact that the Gibbs probability measure exp{—PH) is a signed 
measure because the random coupling is transformed into a complex random magnetic 
field in (17). From fig 1, it is also clear that the mean field solution does not exhibit 
phase transitions at finite temperature. However, there is an essential singularity at 
T = 0, since inserting (24) into (15) and (14) one sees that / ~ exp(l/T) for T ^ 0. 

It is possible to explicitly solve the self-consistency equation when /3 — > since 
(25) becomes 



The zero temperature energy of the mean field solution is Eq = —1.5 while the 
numerical simulations [7] give Eq = —1.404 ± 0.002. In fig 2, we show the free 
energy as a function of T. One sees that entropy is negative at low temperature, 
thus indicating that the solution is unphysical. As a consequence, a better estimate 
of the ground state energy is given by the maximum of /(/?), following a standard 
argument of Toulouse and Vannimenus [8], and one has Eq > max^ /(/?) = —1.468. 

3. BETHE-PEIERLS APPROXIMATION 

The mean field approximation neglects the short range order, that can be taken into 
account by the so-called Bethe-Peierls approximation [9-10] . It is still useful to consider 
the model on the dual lattice and, moreover, it is convenient to work on the internal 



coth(8/?$) = 



(25) 




(26) 
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energy 



= (27) 



instead of the free energy. From (14) and (15) one thus has 



U{P) = -2 coth(2/5) - — ^W(/?) (28) 



with 



W(/3) = ^ (/3/(/3) ) = lim Wn(/3) 
where /(/?) is given by (18) so that the internal energy of n rephcas is: 

{s} i \ a / a 

that in the hmit n — > gives U{P). 

Now comes the key step. Let us sum over all the spin couples but the 4n nearest 
neighbor spins o'2^\ ^3°'\ ^4^^ around the n spin o^q"^ in the numerator of (29). 
In order to simplify the notation we shell indicate in the following all these 5n spins 
as S(5) and the 4n lateral ones as 8(4) 



Noting that the prefactor (^jv)" is a constant depending on n, the expression (29) 
becomes: 



-1 \4n 

w„(« = ((i;^F) E (i n 1 +11'?.'°' nn 

{s(5)} V ^ 1=0 \ a / 1=1 a 



(30) 

where is a function of the 4n lateral spins 8(4). In practice, we are considering 
the 5n free spins 8(5) on replicated crosses which are interacting with the mean field 
generated by the other n{N — 5) ones. 

The ansatz of the Bethe-Peierls approximation consists in assuming that any 
function of the spin configurations such as (30) might be expressed as 

n 4 



^a=l i=l 

10 



where we have introduce the normahzation factor 

z{n,P,fx)= J2 ^n(S(5)) (32) 

{S(5)} 

related to the weight of the S(5) configurations 

= (i)' f 1 + n^r] n f 1 +11-1'"] fi'^'^'""'*'""'"' (33) 

The parameter is a sort of chemical potential representing the energy cost necessary 
to flip the lateral spins in the opposite direction of ao, destroying the short range 
order. In fact, the Bethe-Peierls approximation is also indicated as quasi- chemical 
approximation. 

As a consequence, the internal energy becomes 

Un0)^ -2(5o5^i)n (34) 

Here and in the following indicates the average of an observable A, 

{A)n = z-\n,p,i,) ^ W^n(s(5)) (35) 

{S(5)} 

over the 2^" configurations of the spins on the replicated crosses weighted by W^. The 
chemical potential /x depends on the replica number n and can be determined through 
a self-consistency equation given by the requirements that the average value of the 
dual spin is invariant under translations, i.e. 

(3^o)n = (5i)n ^ = 1, • • • , 4 ; a = 1, • • • , n (36) 

One is interested in the limit n — > 0, as usual. In order to write the self-consistency 
equation in a simpler form, it is convenient to introduce the generating function 

Mh,f^,P) = In J] W^n(s(5))e'^^«^o"' (37) 

{S(5)} 
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so that (36) corresponds to require 



dh 



- 1^ 



(38) 



h=0 



where is the quenched generating function, 



(j){h,^i,(3) = hni — . 

n— s-O n 

The solution of this imphcit equation gives the value of the chemical potential /x*(/3) 
as function of the temperature. The internal energy can then be expressed in terms 
of the generating function as 



U{(3) = - 



2 d(5 



(39) 



In order to obtain the quenched generating function </>, we should perform some 
algebraic manipulations. After performing the sum over the 2'^"'^ configurations {s(4)} 
in (37) we remain with a sum over the configurations Sq = ^q'\ ■ ■ • , o^o"^ 



(a) 



{so} \ a / \ a a J 



E E (fc) 1 + n^^o"' n (2 cosh,*"))'"' n {2.i^r,(-^y e"-''" 



{so} k=0 
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= EEE 

{so} fc=Oj=±l 

where we have introduced the variable 



(a) 



(41) 



for simplifying the notation. Now, the previous sum has been obtained by an annealed 
average over the disorder, i.e. 



EE 2- 

fc=0 j=±l 



k 



A 



(42) 
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where one has 



{so} a 

(2cosh?7)^"'= (2sinhr7)'' e^^° J = (a^j)" (43) 
with rj = /! + (3aQ. Noting that 

1 , 



1 i 1 1 1 — In a"^ = In a 

n— s-O n 

eventually one can write the quenched generating function as 

<^ = ^ E E (^) 1^ E ^ ^i^h'^ ^ (44) 

i=±i fc=o ^ ^ ~ 

Here and in the following we omit to write the constant additive term 4 In 2 in 0. Note 
that such a term disappears in the derivatives. The first sum over j in (44) can be 
performed by a trick. Let us use an auxiliary spin 5q = ±1 with equal probability so 
that 

1 ^ /4\ ~ ~ 

^~2^E(y^)^^ E e''^'^°+'^o) (cosh?7 coshry')'*"'' (sinhry sinhry')'' (45) 

with rj' = fj, + /3?Q. It is easy to realize that if ctq cind 5g have opposite sign, the 
contribution to the sum vanishes. We can limit ourselves to consider the case of equal 
sign, so that (45) becomes 

= ^ E (t) E ^0 cosh^^^-'^) T] sinh^'^ r/ (46) 

Moreover in the limit h ^ 0, one has 

50 6^'*'"° ~ 5o + 2/i (47) 
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It follows that (uo) is 
dh 

with 



° ^—^ cro=±l 



Sk= 3^0 cosh'(^-^) 7/ sinh'S 



(48) 



(49) 



while Ei=i(f^i) is 



djj, 



^ o-o=±l 



with 



Yk{r]) = cosh^^^"'') 7/ sinh^'' r] [2(4 - /c) tanhry + 2k coth 



?7j 



(50) 



(51) 



Inserting (45) and (46) into the self-consistency equation (38), the chemical potential 
jj,* can be obtained as function of the dual inverse temperature Once determined 
the value /i* , the internal energy is given by a derivative of the generating function. 
In particular one has 

1 d(j) 

h=0,ij,*{/3) K=v ~ ' CTo=±l 



d(3 



-Y^tQ'^' E (52) 

fc=0 ^ ^ ~ _- n 



In fig 3 we show fi* / (3 as a function of /?, where one observes that for /? — > oo (limit of 
high temperature T of the original lattice), fi* / P 3. This can be understood noting 
that each one of the lateral free spins of the cross interacts with three other spins so 
that, at zero dual temperature the energy lost in a flip is exactly equal to 3. Fig 
4 illustrates the graphical solution of (38) by plotting 



lim (^{ao)n - \ XJ(^z)n j 



(53) 



as function of \x at three different temperatures. The solution \x* of the self-consistency 
equation is given by the intersection of the function with the horizontal axes. We look 
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only for real solutions. At large there exists only one solution. However, fig 4 
shows that at low j3 two solutions appear, and for (3 < 0.031 there is no real positive 
solution. The internal energy U given by the Bethe-Peierls approximation together 
with the annealed energy Ua = — 2tanh(/3) and the mean field result are plotted as 
function of the temperature in fig 5. One sees that for /3 < 0.05, i.e. T < 0.667.., 
the energy increases at lowering the temperature while the chemical potential //* 
decreases, indicating that the Bethe-Peierls solution becomes unphysical. The ground 
state energy can be estimated by the minimum value assumed by the internal energy, 
i.e. 



It is extremely close to the numerical estimate of [7] Eq = —1.404 ± 0.002. 

Finally, we want to mention that it remains open the problem to understand 
whether, with an appropriate ansatz of replica symmetry breaking, one can obtain a 
solution of the self-consistency equation for /?< 0.03. It is indeed well known that 
mean field approach can give phase transitions, even when they are absent in the 
original model. 

4. IMPROVED BETHE-PEIERLS APPROXIMATION 

In order to improve the Bethe-Peierls approximation in the framework of the replica 
method, we introduce a second variational parameter 7, beyond the chemical potential 
fi, that puts in interaction the different replicas. In other terms, we replace the 
(standard) ansatz (31) with 



Eq rmnU{p) = U{(3 = 0.05..) = -1.3975 

/3 



g{n,/3,n,j) 



(54) 



a i=l 



where we have introduced the normalization factor 




(55) 



{S(5)} 
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related to the weight of the 8(5) configurations 



^n(S(5)) = WnC 



(«)~(/3) 



a: 'a 

a>f3 I 



(56) 



with Wn given by (33). We shall indicates the average of an observable A over this 
new normalized weight by << A »n- 

The two variational parameters are determined by the couple self-consistency 
equations obtained in the limit n — > by 



<< (To »n = « CTi »n 



(57a) 



(576) 



with i = 1, • • • , 4 and a, P = 1, ■ ■ ■ ,n. 

The ansatz proposed here is related to a hypothesis of existence of a glassy phase. 
In fact, one can apply this approximation scheme to the the solution of the random 
coupling Ising model directly on the original lattice in d-dimensions [11]. In this 
case, after performing the limit d — > 00, one obtains the Parisi solution [2] of the 
Sherrington-Kirkpatrick model. 

Following the same idea of the previous section, let us introduce the generating 
function 

(58) 



V',(/i,£,^,7,/?) = ln J2 ^'n(s(5))e'^^«"°" +^E„>,-i"- 

{S(5)} 

SO that (57a) and (57b) correspond to require 



and 



dip 
dh 



dip 
'di 



- 1^ 
- 



(59a) 



h=Ql=Q 



(596) 



h=Ql=Q 
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where is the quenched generating function, 

V'(/i,//,7,^)= hm^. (60) 

n— >o n 

The solution of this imphcit equation gives the values and 7* as function of the 

temperature. The internal energy can then be expressed in terms of the generating 
function as 

(61) 



^ ' 2 d(3 



/i=0,^=0,M*(^),7*(^) 

Let us now simplify as much as possible the self-consistency equations. It is 



convenient to use the standard gaussian identity 



In 



in order to write the generating function as 



(62) 



(47 + £)n , / 5 fJjdxi 

{so} \ a / \ cx a / J 

with 

^(a)=^(a)+^.^ (64) 

Using the same 'algebraic' strategy that in the previous section leads to (44), we obtain 
the quenched generating function as 



4 



£ 1 /" TT dXi 

»=0 " 3i,3i,33,]A=±l 



ln[ J] 5oe^('^+^^°^''° J|(cosha;i)i+^'* (sinha;^)^-^'* I (65) 
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with 

LOi = r] + Xi^/y ^ (3ao + fi + xi^/^. 

The constant additive term 4 In 2 is again omitted. 

By derivating '0 one has the self-consistency equations (57) for and 7* in terms 
of the sum of 5 gaussian integrals. A careful analytic and numerical study of these 
equations might enlighten the nature of spin glasses in low dimensional systems. 

5. CONCLUSIONS 

Let us briefly summarize our main results: 

(1) formulation of the random coupling Ising model on the dual lattice made of square 
plaquettes. The dual model has signed Gibbs probability measure as the random 
coupling is transformed into a random complex magnetic field. 

(2) Application of the mean field approximation to the two dimensional dual model 
in the framework of the replica method. We find the solution using a replica 
symmetry ansatz, obtaining a good estimate of the ground state energy: Eq = 
—1.468, compared with the numerical result of [7] Eq — —1.404 ± 0.002. 

(3) Application of the Bethe-Peierls approximation. It gives a very accurate estimate 
of the ground state energy (Eq = —1.3975) although it becomes unphysical below 
P = 0.05 and there is no real solution of the self-consistency equation for P < 0.03. 

(4) Improvement scheme of the Bethe-Peierls approximation by considering a second 
variational parameter that puts in interaction different replicas of the dual model. 
There are still many problems that remain open in our approach. As major issue, 

it would be interesting to understand whether the ansatz proposed in section 4, or 
other similar assumptions lead to a solution of the self-consistency equations of the 
Bethe-Peierls approximation at low (3, that allow one to decide whether a transition 
to a glassy phase is present at low dimension. 
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The dual transformation is indeed a very powerful tool for determining the critical 
temperature in non- disordered systems and our method might give some results in 
this direction. Least but not last, we plan to find a cluster expansion scheme that 
permits to improve the mean field in a systematic way. 
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Figure Captions 

Fig. 1 Graphical solution of the implicit equation (25) at T = /3~^ = 1 corresponding to 
P = 0.136... The full lines are coth(8/3$) versus $ and the straight line $ = $. 

Fig. 2 Annealed free energy given by (4) (dashed line) and the mean field solution 
(full line) versus temperature T — I3~^ . The dashed lines are the Maxwell 
constructions obtained by imposing that the free energy is a monotonous non- 
decreasing function of T. The annealed solution estimates a ground state energy 
Eq > —1.559; the mean field solution gives Eq > —1.468; the numerical result of 
[5] is Eo = -1.404 ±0.002 

Fig. 3 Chemical potential as function of the dual inverse temperature f3. 

Fig. 4 The Graphical solution of the self-consistency equation (36) is given by the 
intersection of (53) with the horizontal axes. We illustrate three different cases: 
P = 0.1 (fuU line), /3 = 0.05 (dashed line) and /3 = 0.028 (dotted line). 

Fig. 5 Internal energy given by the annealed approximation, i.e. Ua = — 2tanh/3 (dashed 
line), by the mean field approximation (dash-dotted line) and by the Bethe-Peierls 
solution (full line) U (T) versus temperature T = (5~^ . The dotted line is obtained 
by imposing that the Bethe-Peierls internal energy is a monotonous non-increasing 
function of T. 
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